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Abstract. In this paper, we have made an accurate investigation of proton acceleration in GRBs and we have 
predicted a possible signature of cosmic rays, in a sufficiently baryon-loaded fireball, via GeV 7-ray emission 
produced by 7r°-meson decay. If two ungrounded assumptions are removed, namely, Bohm's scaling and a slow 
magnetic field decrease, the usual Fermi processes are unable to generate ultra high energy cosmic rays (UHECRs) 
in GRBs. We propose to develop another scenario of relativistic Fermi acceleration in the internal shock stage. 
We present the results of a realistic Monte-Carlo simulation of a multi-front acceleration which clearly shows the 
possible generation of UHECR. The amount of energy converted into UHECRs turns out to be a sizeable fraction 
of the magnetic energy. 

Key words, gamma rays: bursts - ISMxosmic rays - acceleration of particles 



X 



1. Introduction 

Gamma Ray Bursts (hereafter GRBs) have been con- 
sidered as promising sources of Ultra High Energy Cosmic 
Rays (UHEC Rs) either t hrough the external ultrarela- 
tivistic shock l | Vietril ll995) or through the internal shocks 
l)Waxmanlll995|h The acceleration at the external shock 
does not achieve this goal when the external med ium 
is a s tandard interstellar medium ijGallant fc Achterberd 
119991) . As for the internal shocks, first and second order 
Fermi processes would achieve this goal with the extreme 
assumption of a cosmic ray mean free path of the or- 
der of its Larmo r radius (so-called Boh m's scaling). In 
a previous paper ilGialis fc Pelletieife oO.Ti , we have shown 
that such scaling is physically inconsistent for two rea- 
sons; first, the law is no t validated by numerical simula- 
tions ()Casse et al.ll200l|) and second, even if inappropri- 
ately used, it would lead to tremendously fast magnetic 
energy depletion in the fireball. We also have considered 
a more realistic Kolmogorov scaling but a severe expan- 
sion loss limitation does not make it possible. These ar- 
guments are developed in Sects. [3 and [21 of this paper, 
together with the important issue of the magnetic field 
profile. Because both the global energy budget of the fire- 
ball and the Hillas criterium are in favour of the genera- 
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tion of UHECRs, we have made a deta iled investigation of 
a diffe rent Fermi pro c ess pr oposed bv lPelletier fc Kersalel 
(2000) and IPelletierl lIlQQQh . In this kind of Fermi pro- 
cess, reviewed in Sect.0J the sheets invoked in the internal 
shock model are considered as relativistic hydromagnetic 
fronts that scatter cosmic rays. Thus, cosmic rays undergo 
a relativistic Fermi acceleration, similar to a kind of sec- 
ond order process through multiple interactions with all 
the fronts in the relativistic wind. Actually, this is not a 
second order Fermi process because no expansion is al- 
lowed in the relativistic regime, the energy jump being 
large at each scattering. The generation of cosmic rays 
resulting from these multiple scatterings is very different 
from generation obtained by a superposition of the indi- 
vidual first order contributions of each shock. 
In Sect.|31we describe a realistic Monte-Carlo simulation of 
this process, taking into account several dynamical param- 
eters in a conical expanding fireball, for several intensities 
and profiles of the average magnetic field. The results are 
presented in Sect. El they suggest a good efficiency of the 
process of acceleration of UHECRs and give a prediction 
of the time at which they are created during the internal 
shock stage. Finally, we establish a possible diagnosis of 
cosmic rays in GRBs through 7-ray emission by 7r°-meson 
decay which would be easily detected by HESS observa- 
tory, 5@5 experiment and GLAST. 
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2. Preliminary considerations 



2.2. Internal shocks 



The goal of the generation oif UHECRs in GRBs is to 
reach a maximum energy close to the confinement energy 
limit above which the Larmor radius becomes larger than 
the correlation length of the magnetic field, namely 



ZeTBi 



(1) 



In Eq. Q, B is measured in the co- moving frame; the 
correlation length, £ c , also measured in the co-moving 
frame, stretches with the expansion of the fireball. This 
length keeps the same order of magnitude as the char- 
acteristic size of the subshells and follows the same scal- 
ing, where the Lorentz factor is related to the velocity (3q 
through T = (31 This is the reason why we as- 

sume l c ~ AR ~ r/r taking r as the distance of a fluid 
element from the explosion center in the observer frame; 
r is the bulk Lorentz factor after the saturation radius. 
Hereafter, we suppose that T reaches the maximum value 
r\ = E/M b c 2 where E is the total energy emitted and M b 
is the total baryonic mass ejected. Then, Equ. Q becomes 



e c i c~ 10 2 



\ 10 6 G J 



io- 



J pc 



eV. 



(2) 



The Fermi acceleration by internal shocks starts at a dis- 
tance of the order of r b = t^tq (about 10 12 cm) which 
is currently defined as the broadening radius i.e the dis- 
tance after which the shells broaden. Thus, in this model, 
a magnetic field of at least 10 6 G is necessary at this ra- 
dius to achieve UHECR production. Since the product B r 
decreases like r x ~ a (with 1 < a < 2), super-GZK protons 
could be produced between r b and 100 r b for a = 3/2 and 
between r b and 10 r b for a = 2. We note that no genera- 
tion of UHECRs is possible for B smaller than a few 10 5 
G in the internal shock model. Regarding the Fermi accel- 
eration at the external ultrarelativistic shock, its intensity 
must reach 100 G at the deceleration radius (~ 10~ 2 pc) 
corresponding to the equipartition value with the pressure 
of the shocked plasma. Let us examine this case first. 



2.1. External shock 

The Fermi acceleration of supra-thermal parti cles at 
the ex ternal ultrarelativistic shock was sugges ted bvlVietri 
1995) and it has recently been reviewed bv IVietri et al " 



2003). As mentioned before, a 100 G magnetic field 



is required to obtain UHECRs. If the usual interstellar 
medium is considered, with a /iG magnetic field, even with 
the shortening of the acceleration time due to the ultra- 
relativistic shockfront, namely t u ~ t^/Ts, as shown by 
Gallant & Achterberg (1999), UHECRs cannot be pro- 
duced. Moreover the amplification of the pre-existing cos- 
mic ray distribution does not provide enough UHECR 
flux. The only hope is through the consideration of a differ- 
ent external medium, as for instance a pulsar wind b ubble 
such as considered bv lGallant fc Achterberd l|l999|) . 



The internal s hock model has been designed 
l|Meszaros et al.lll993|) to account for variability observed 
in the light curve. Indeed millisecond variability, as ob- 
served in some GRB light curves, cannot be explained in 
the framework of a single sheet colliding with the interstel- 
lar medium. The most interesting feature of the multiple 
sheet model is that the time variability at the origin of the 
flow is observable thanks to the following scheme. Let us 
consider two sheets that leave the central engine separated 
by a time interval At, respectively with Lorentz factor Ti 
and T2- A collision occurs when T2 > T% at a timr that 
can be expressed by 



2r 2 r 2 



r 2 
1 2 



At. 



Assuming an instantaneous shock pulse, some time 
spreading, At s — t c /2T 2 , is observed and thus this 
duration corresponds to the original time interval. The 
shortest variabilities (~ 1 ms) correspond to the typical 
time scale associated with the size of a black hole of a few 
tens of solar masses (namely Tq/c). Such a first collision 
takes place around the distance r b . Longer variations 
correspond to collisions at a more remote distance up to 
a maximum distance determined by the duration of the 



flow At w . This maximum distance is 



r b cAt w /r , 



with cAt w /ro ~ 3 x 10 3 (At w /ls), which gives the proper 
length of the flow in the co-moving frame £0 — (3cT At w . 
The duration of the flow during the internal shock 
phase is therefore At max ~ (rb/ro)At w ~ rj At w . 
This phenomenon is observed during a time in- 
terval shortened by the propagation effect, namely 
At obs = (1 - [i)At max ~ At max /2r] 2 ~ At w . The previ- 
ous value of r max determines the range of the acceleration 



process that extends from r& to 



which is not far 



from the deceleration radius r& ~ 10 cm. 

The Fermi acceleration (first or secon d order) in the 
internal shock model is usually considered \ Waxmanl 1 995|) 
as mildly or sub-relativistic with a characteristic time pro- 
portional to the Larmor radius (Bohm sc aling). However, 
in a recent paper (|Oialis fc Pelletierir2 003) , we have shown 
that this assumption is not realistic regarding the mag- 
netic energy depletion time. Moreover, the Fermi acceler- 
ation time depends on the mean free path, £, of the par- 
ticle in an irregular magnetic field. This length depends 
on two other lengths, namely the Larmor radius, tl, and 
the correlation length, £ c : for a turbulence spectrum of 
magnetic perturbations that follows a power law of index 
(3 the following law, which is known in weak turbulence 
theory, has been extended into the regime o f strong tur- 
bulen ce and large rigidities such that < £ c 1 Casse et alJ 
12001: 



Vt U 



i-0 



(3) 



D. Gialis and G. Pelletier: Which acceleration process for UHECRs in GRBs ? 



3 



where rj t = ^rg^ . The Bohm scaling i ~ tl, which holds 
for electrostatic turbulence, does not apply for purely 
magnetic irregularities on large scales; no theory nor nu- 
merical simulation has confirmed Bohm's conjecture. The 
Bohm estimate corresponds only to the specific case where 
the magnetic field is totally disorganised and the Larmor 
radius is as large as the correlation length, which is not 
the case in GRBs. 

Accord ing to a Kolmogorov sca ling with (3 — 5/3, we have 
shown (|Gialis fc Pelletierl l2003) that GRBs are unable to 
produce UHECRs with this acceleration process because 
of a strong expansion limitation in energy 




GeV, (4) 



where kq is the ratio of the acceleration time over the 
Larmor time for a Larmor radius that equals the corre- 
lation length of the magnetic field. This limitation, mea- 
sured in the co-moving frame, is more severe than the 
synchrotron one and suggests that we have to consider 
another type of process to achieve high energy. This is the 
purpose of Sect. [J] 

2.3. The magnetic field evolution 

The flow that reaches the beginning of the acceleration 
range is disconnected from the source if cAt w < r\>. In that 
case, the sheets that compose the flow start to expand 
with car rying a magnetic field t hat decays like r -2 , as ar- 
gued by iMeszaros et alJ l|l99.3fo . This is the most severe 
assumption on the evolution of the magnetic fi eld in the 
attem pt to accelerate particles. Some authors (jWaxmanl 
120031) claimed that it decays like r _1 , because this would 
be consistent with the synchrotron interpretation of the 
GRB spectrum; this would require an amplification of the 
magnetic field in the vicinity of the shocks. We think that 
this would occur only in narrow layers behind the internal 
shocks and would be beneficial for electron acceleration. 
As for the high energy protons, especially those candidates 
for UHECRs, they necessarily explore a major part of the 
flow because of their large Larmor radii, and they experi- 
ence a magnetic field that decays faster than r _1 , unless 
the flow is still connected with the source and carries a 
current associated with a toroidal field. In computations, 
for the sake of compromise, we have considered a field that 
decays like r~ 3 / 2 , which means that the magnetic energy 
is conserved in the shells. However, to demonstrate the 
capability of GRBs to produce UHECRs, we will also in- 
vestigate the conservative case of a magnetic field decaying 
like r~ 2 . 

Because we will need some magnetic field intensity at 
the beginning of the acceleration range, we need to assume 
Bo ~ 10 15 G at the black hole scale r$ ~ 10 7 cm. However, 
for the work presented in this paper, we do not consider 
a situation where the magnetic field would dominate dur- 
ing the adiabatic expansion. The case of a magnetic field 
dominated expansion is interesting from the point of view 



of particle acceleration, but requires different modeling, 
postponed for future work. 

3. The distribution function resulting from 
internal shocks 

In this paper, we will introduce an acceleration pro- 
cess additionnal to the usual internal shock acceleration. 
However we will consider this additionnal process as 
operating above 1 GeV and below the expansion cut 
off mentioned in Sect. 12.21 (see Eq. Q). This yields an 
e~ 2 -spectrum in the frame TZt moving with a Lorentz 
factor r, in which the distribution function is almost 
isotropized. In fact, we will see later on that the highest 
energy part, especially that part accelerated by the 
additional process, is not fully isotropized. Because we 
will present numerical results in the observer frame (or, 
more properly, in the stationary frame of the source TZq), 
it is useful to show how this truncated distribution looks 
in this frame. 

In the co-moving frame IZt , the distribution is almost 
isotropized by turbulence, at least for the lower energy 
part, and the phase-space distribution function can be 
written 

r{tco^ co )=Ke-J (5) 

with K — 47T ^y f° r 1 < e co < 7m and it is null 
otherwise. The Lorentz factor 7 m corresponds to the 
cut-off energy due to expansion. 

In the stationary frame IZo, the distribution is highly 
anisotropic, since the proton energy, e, is such that e = 
5(p)e co where S(fi) = 1/ (T (1 — po fi) ) . The angular distri- 
bution is such that 

^ = 2L 2 (1 - p^) 2 ' (6) 
where L ^> 1 implies values very close to 1 for /1. 

The Lorentz invariance of the phase-space distribution 
function leads to an energy distribution function / such 
that 

r 1 IttKf- 2 r 2T 

f(e)^2ne 2 /(e, M )rf M = / P(S) 6 2 dd(7) 

J -i Pol J1/2T 

with P{8) = 1 for e/j m < S < e and otherwise. 

In the frame TZq, one can summarize the effect of the 
anisotropy on the shape of the energy distribution function 
produced by the internal shocks as follows (see appendix) 

f e for e < 2L 

/"(e) cx <^ e~ 2 for2r<e<2L 7m . (8) 
[ otherwise 

As we will see in Sect. |SJ a numerical simulation gives 
an illustration of the previous results about this distribu- 
tion function resulting from a first stage of internal shock 
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acceleration. It constitutes the "initial" condition of an 
additional Fermi acceleration process that we consider in 
the next section. 

4. Relativistic Fermi acceleration by multiple 
fronts 

In this section, we present the properties of hydromag- 
netic fronts in GRBs and we show that the Fermi accel- 
eration process resulting from scattering off relativistic 
magnetized fronts constitutes an interesting solution for 
UHECR generation. 

4.1. Characterization of the hydromagnetic fronts in 
GRBs 

The internal shock model accounts for the light curves 
of GRBs that display fast variations, sometimes on a 
millisecond scale. This typically requires a number N s 
of sheets around 10-100, with Lorentz factors Tk (k = 
1, 2, Ns) between 10 2 -10 3 w ith respect to the frame Ho 
l|Daigne k Mochkovitcblll998|) . This means that they are 
mildly relativistic in the co-moving frame TZr with Lorentz 
factors 7fe ~ 2, and we will refer to an average value 7*. 
For instance, assume that half of the sheets run forwards 
at a velocity /?* with respect to TZt, and half of them 
run backwards at a velocity — then we have a lower 
Lorentz factor T~ = and an upper Lorentz factor 

r+ = 27*r with respect to IZo- 

The cosmic rays can be scattered by these sheets 
provided that they carry a sufficiently high magnetic 
field. Indeed they are scattered by magnetic sheets if 
and only if their Larmor radius is smaller than the 
thickness Sk of the magnetic profile of the considered 
sheet. Particles of larger Larmor radius cross the front 
with just a small pitch angle jump of order Sk /tl. We 
assume that Sk is of the order of the thickness of the 
sheet and follows its evolution. Therefore, we assume an 
average thickness 5 ~ r/T 2 measured in the stationary 
frame TZq. Consequently, particles of energy larger than 
Z e r B+ 5, where B+ is the average magnetic field in a 
sheet located at radius ?*, have a negligible interaction 
with the magnetic front. Thus, it is important to specify 
the interaction limit, e^, for each sheet along the flow. We 
adopted the following scheme : we assume a prescribed 
value of the magnetic field .B+(r&) for each sheet at the 
distance rj, and then any sheet carries a magnetic field 
of intensity decreasing with the distance r like r~° with 
a £ [1,2] (we will look at the case a = 3/2 and a = 2 for 
the reason explained in Sect. I2.3|l . Now in the frame of 
each sheet, the particles are scattered with no preferred 
direction and we can adopt either a mirror reflection or a 
random scattering. Both will give similar energy gains as 
will be seen later on; however we will choose the random 
scattering for the numerical simulation. 
Very high intensities of the magnetic field are allowed 
if we consider that the fireball results from the merg- 
ing/collapse of a compact neutron star, namely 10 15 



G. For instance, a decrease of the magnetic field with 
a = 3/2 and an average Lorentz factor of 300, lead to 
J3*(r & ) ~ 3 x 10 6 G for 100 sheets. We emphasize that 
this is a maximum value for a given number of sheets. 
Thus the specification of the magnetic field is crucial, 
for it controls both the acceleration through scattering 
of particles having an energy lower than the interaction 
limit £k oc r 1_Q and the spreading and decimation of the 
scattering sheets, responsible of both acceleration and 
expansion loss. 

Because of the forward and backward, mildly relativis- 
tic, motions of the sheets with respect to the co-moving 
frame IZr, no isotropic distribution function can be built, 
only an axisymmetric and mirror-symmetric distribution 
is obtained. Assuming for instance a random scattering 
scheme (i.e. tendency of isotropization by each sheet), we 
obtain an angular distribution in this frame which is close 
to a superposition of two anisotropic distributions gen- 
erated through the isotropic scattering by both forward 
(moving at +/?*) and backward (moving at — /?*) pertur- 
bations : 

1 1 + a 2 y 



p(m) 



272 (1 - pi^y 



(9) 



4.2. The relativistic Fermi process 



The efficiency of the Fermi process with relativistic 
scattering fronts for the ac celeration of sup ra-thermal par- 
ticles was emphasized by IPelletierl l)l999fl . This process 
consists in an interaction of a relativistic particle having 
an initial 4-vector momentum (ei,pi) in the stationary 
frame TZq with a magnetized front of velocity (3k (with 
corresponding Tk). We denote by primes physical quanti- 
ties in the magnetized front frame IZk and we define the 
angle between (3k and p\ by its cosine \i\ just before the 
interaction. The Lorentz transformation gives us the fol- 
lowing equations : 



p[ ~ T fe (1 - fikHijpi 



1 

Mi 



fj-i - (3k 



1 - (3kH\ 

and assuming e% ~ p\ c (> 1 GeV). 



(10) 
(11) 



The scattering does not change the particle energy in 
the magnetized front frame and the pitch angle (cosine) 

is randomly changed into fi' 2 . Denoting by 2 the phys- 
ical quantities after the interaction, the reversed Lorentz 
transformation leads to 



P2 



M2 



Y 2 k (l+(3k^){l-(3 k ^)pi 



/4 



(3k 



1 + M 



(12) 
(13) 



However to estimate the energy gain, it is necessary to 
calculate this in the frame where the distribution function 
is even in terms of the cosine of the pitch angle, namely in 
7£r- The following unprimed quantities are now measured 
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in TZ-p. Therefore the energy gain per scattering spreads 
from (1 — /?&)/( 1 + At) — 1/47| (for rear-on collision) 
to (1 + /3fc)/(l — Pk) — (for head-on collision), the /3 
being the velocities measured in IZp . Since the distribution 
function is even in [i in 7?-r, the average energy gain per 
interaction is ~ 7 2 for both the random scattering and 
the mirror reflection schemes, 7* being the average value 
of the 7/c. For example, a proton which initially has an 
energy of 100 GeV could achieve an energy of 10 10 GeV 
after about 12 interactions with 7* = 2, whereas 6 Fermi 
cycles would amplify its energy by a factor 64 since the 
gain per Fermi cycle remains smal ler than 2 (after the first 
one) l|Gallant fc Achterberdll99^) . Therefore the interest 
of multi-front interactions compared to the acceleration in 
a single relativistic shock is that the full energy gain ~ 7 2 
can be obtained at each scattering with crossing fronts, 
and could generate UHECRs if the cosmic rays undergo 
enough scatterings off the magnetized fronts. We will see 
further on that this required number of scatterings puts a 
constraint on the number of sheets. 

5. Description of the numerical simulation 

5.1. Initial conditions 

We have considered a conical box representing the 
ultrarelativistic flow around the broadening radius rt, 
(~ 10 12 cm). In the stationary frame TZq, the initial 
box height, corresponding to the flow length, is given 
by cAt w and the transversal size of the flow is equal 
to y/n/irr with fl ~ An/ 500, the opening angle of the 
flow. Using a Monte-Carlo simulation, we randomly put 
in the flow a set of identical (same mass) baryonic lay- 
ers (or sheets) with an uniform Lorentz factor distribu- 
tion between 10 2 and 10 3 considered by s everal authors 
l|Daigne fc Mochkovitchlll99St IPirarJ l200fl to reproduce 
the GRB light cur ves from internal shocks . Following pre- 
vious simulations ^Kobavashi et al.lll997j) . we consider a 
number of layers of about several times ten. For a layer 
with a Lorentz factor of Tk at a distance r the width is 
equal to r jY\ and the transverse radius is the same as 
that of the flow. 

The isotropic suprathermal proton population is initially 
injected in each layer with an energy spectrum in the co- 
moving frame proportional to e' 2 between 1 GeV and a 
cut-off energy 

^'° 4 (w) GeF ' < 14 > 

according to the strong expansion limitation in energy (see 
Eq. Q). In the stationary frame TZq, a Lorentz transfor- 
mation provides the energy, e, of a proton 

e = r k (l±[3 k p c r )e' , (15) 

where is the radial velocity of the proton in the co- 
moving frame of the sheet and d its energy. One can note 
that the energy e' of a proton in the co-moving frame IZk 
can be amplified by a factor 2T k in the stationary frame. 



For each proton, we initially define a random pitch angle 
cosine as defined in Sect. |3] (see Eq. JBJ). We have plot- 
ted in Fig. □ for B*(r b ) = 10 5 G and B*(r 6 ) = 10 6 G, 
the initial energy repartition of the proton population, in 
the observer frame, which depends on three independent 
random parameters namely T kl f3£ and e'. 

5.2. Temporal evolution 

We describe here the temporal evolution in our Monte- 
Carlo simulation. The Lorentz factor distribution leads 
to collisions of different sheets during the numerical pro- 
cess between n, and 300 (about for r c i = 10 15 cm 
or rd/30 for — 10 16 cm) which is sufficient to deter- 
mine the fraction of UHECRs according to the confine- 
ment limit (see Sect. [2J - We have considered the simple 
case where the result of a collision of two identical sheets 
is a single sheet having a Lorentz factor defined by \JT\ T2 
where Ti and T2 are the Lorentz factors of the sheets be- 
fore colliding. The total number of sheets decreases with 
time and this reduces the interaction probability with pro- 
tons in the flow even if the sheet width and the transversal 
radius increase like r in our model. 

The first collision of a sheet defines the time from which its 
protons are released into the flow and can interact with all 
sheets. After a proton-sheet interaction, which is spatially 
defined by comparing respective positions, we randomly 
assign a new pitch angle for the proton and its energy is 
changed according to the equation (15} . We assume that 
all interactions are instantaneous. Indeed, the magnetic 
scattering duration inside a sheet, r s , can be neglected; 
for a proton with energy e' in the co-moving frame, it can 
be expressed by 

r.-lST.xlO-*^)" 1 ^)., (16) 

and it is very small compared to the mean free path time 
r/cN s > 1 s. 

The proton energy limitation for an interaction (or mag- 
netic scattering) is given by the interaction limit calcu- 
lated in each layer (see Sect. |21 and paragraph 14. 1(1 : for a 
given proton, the higher its energy is, the more transpar- 
ent the medium is to it and its escape probability becomes 
important. This limitation is far more important than the 
transverse escape. We also have neglected the synchrotron 
losses for protons because the synchrotron limitation, in 
the co- moving frame IZr, is over 10 5 GeV at and next, 
increases like r 2 for B+ oc r~ 2 and r 11 / 8 for B* oc r~ 3 / 2 
llGialis fc Pelletierll2003h . 

6. Numerical results and discussion 

6.1. Generation of UHECRs 

The numerical Monte-Carlo simulation clearly shows 
that a sizeable part of cosmic rays is produced by the 
relativistic Fermi process stretching the initial distribution 
tail up into the UHE range (see Fig. [2J) ; for any magnetic 
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Log(Energy) (GeV) 



Fig. 1. The proton repartition in the stationary frame TZq calculated with the simulation and before the additionnal 
acceleration Fermi process : this proton population is just the result of a statistical superposition of the contribution 
of every sheet. Top for a magnetic field at rb of 10 5 G. Bottom: B±(rb) = 10 6 G. According to the initial distribution 
function calculated in Sect.|3|and the cut-off energy CEq. (| 1411 1. the higher the magnetic field, the more the distribution 
spreads over high energy. 



field intensity higher than 10 4 G at rb, the relativistic 
Fermi acceleration by multiple fronts leads to a production 
of UHECRs with a fraction larger than 1CP 6 for a given 
initial proton distribution as prescribed in Sect. [5] Table 
n summarizes the results for different values of B+(rb) and 
of the flow duration. We have chosen a typical number 
of sheets, N s , of 50 at the beginning of the process and, 
first, a magnetic field that decreases like r -3 / 2 , second, a 
magnetic field that decreases like r~ 2 . 

The determination of the time of UHECR creation 
which has been made for different values of the wind dura- 
tion and the magnetic field (see Fig.0} shows that UHECR 
generation mainly happens at the very beginning of the 
internal shock phase; more precisely, UHECR generation 
can be achieved before ~ 300 rj (or ~ 10 4 s) which corre- 
sponds to only a few percent of the internal shock phase 
duration (> 10 5 s). There are two conditions that gov- 
ern the ending of the multi-front acceleration: first, the 
maximum value of the interaction limit assigned at the 
radius rb through the value of the magnetic field there; 
and second, the decrease of the interaction energy limit 
with distance. Indeed, the lowest interaction energy limit 
(see Sect. EJ to get UHECRs is about 10 19 eV at r b for 
B±(rb) — 10 4 G; the difference between this energy and 
the UHE range (> 10 20 eV) is almost at the maximum 
possible value of the energy jump (the maximum gain be- 
ing about 16), and then the gap quickly becomes too high 
with time to be jumped by a proton in a single proton- 



sheet interaction, since the interaction limit decreases like 
efc cx r x ~ a . However, we have shown (see Sect, 0} that 
only a small number of interactions (> 10) is necessary 
to achieve the UHE range. The simulation directly pro- 
vides this number per proton: whatever the flow duration 
between 0.5 and 10 s and for a magnetic field of 10 4 G 
at rb, the results confirm that a sizeable fraction of pro- 
tons (> 2.0 x 10~ 3 ) undergoes more than 10 interactions 
with the hydromagnetic fronts and, consequently, can eas- 
ily achieve the UHE range (see Fig. 0J. Moreover, Tabled 
indicates that the energy amount converted into UHECRs 
at the beginning of the Fermi process is high enough to 
constitute a sizeable fraction of the total magnetic energy. 
Finally, a magnetic field decreasing in r~ 2 leads to a lower 
resulting fraction of UHECRs but a sizeable one again: for 
instance, for a magnetic field 5* (?-(,) = 10 4 G, this fraction 
is about 8.0 x 10~ 6 with At w = 10 s and 2.6 x 10~ 5 with 
At w = 0.5 s. 

Let us now give arguments concerning the number of 
sheets required to produce UHECRs. The numerical sim- 
ulation shows that the acceleration stops because of the 
decimation of the sheets due to expansion (the mean free 
path increases like l s cx r) and mostly due to decrease 
of the interaction energy limits (e& cx r _1 for B+ cx r~ 2 ) 
which rarefies the scattering medium for increasing en- 
ergies. Since £ s /r is given by the number of sheets N s , 
namely £ s /r — 1/N S , a simple estimate of the minimum 
number of sheets can be given. At each scattering, a parti- 
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Fig. 2. The initial proton distribution (left figure) is strongly modified by the Fermi acceleration process in the high 
energy range where an UHECR population appears (right figure) with B*(rb) = fO 5 G, At w = 10 s, N s = 50. 



Flow duration (in s) 


Ma 

10 4 G 


gnetic field a 
10 5 G 


10 6 G 


0.5 


5.4 x 10~ b 


1.5 x 10~ 4 


3.6 x 10" 4 


1.0 


6.3 x 10~ b 


1.9 x 10" 4 


4.1 x 10" 4 


5.0 


5.1 x 10~ b 


1.8 x 10~ 4 


4.8 x 10" 4 


10.0 


2.0 x 10~ b 


1.2 x 10" 4 


3.2 x 10~ 4 



Table 1. The UHECR (> 10 10 GeV) fraction generated by relativistic Fermi acceleration in the Monte-Carlo simu- 
lation. A number of 50 layers is considered here and the flow duration varies from 0.5 s (short GRB) to 10.0 s (long 
GRB). The magnetic field decreases as r~ 3 / 2 . 



cle migrates towards a region of lower interaction limit and 
roughly, after n scatterings e s (n) = (1 — q) n e s (0), where 
q = (a — 1)/N S and e s denotes the local scattering energy 
limit, rather than the interaction limit associated to a spe- 
cific sheet. So the acceleration stops when the number of 
scatterings is such that 7 2 ™eo ~ (1 — q) n e s (0), which leads 
to the following number of scatterings: 




Since we expect that the energy reached after n c scatter- 
ings is not far from the maximum value of the interaction 
limit, it requires n c q < 1, and thus 




To get particles that will reach an energy at least equal to 
the interaction limit at with 7* = 2 and B+(rb) = 10 4 
G, we need more than 13 to 17 sheets. 

Since the process tends to produce cosmic rays that 
reach an energy close to the local scattering limit, the in- 
tegrated distribution function reflects the evolution of the 
scattering limit, which is directly given by the magnetic 
field decay. For an intensity of the magnetic field decreas- 
ing like r -2 , the scattering limit decreases like r _1 , and 
for a continuous flow during At w , the cosmic ray density 
decreases like r~ 2 . Therefore the integrated energy distri- 



bution function is 

-^T- = J <5(e-e 5 (r))n^r2r 2 dr cx J S(e-e s (r)) cx e~ 2 . 

(19) 

Thus we obtain an e 2 -spectrum over 4 decades, because 
the range of scattering limit energy stems from the accel- 
eration range that extends from rj to ~ r^. 

6.2. Diagnosis of Cosmic Rays in GRBs 

6.2.1. The two types of neutrino emission 

There are two possible types of neutrino emission du e 
to the generation of cosmic rays llGialis fc Pelletierl feoOS) . 
The most expected one results from the p7-process where 
the total energy radiated by neutrinos in the co-moving 
flow is given by 

El~ao£Mog(^) e;, (20) 

V ith J 

where the threshold is at jth ~ 10 4 - It can be easily 
established that even this improved acceleration process 
does not provide a sufficient increase of the neutrino 
emission. Taking account of the collimation, the number 
of p7-neutrinos that cross a km 2 -detector, coming from a 
GRB at lOOMpc, is about 10 4 . 

The other possibility is low energy emission between 
5 GeV and 150 GeV resulting from pp-collisions. We pre- 
dicted significant emission due to the standard Fermi ac- 
celeration at internal shocks when rj is smaller than rj±/2 
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Fig. 3. The proton repartition as a function of the number of interactions (or scatterings) for a magnetic field at of 
10 4 G and for a wind duration At w = 10 s: a sizeable fraction of protons (> 2.0 x 10~ 3 ) undergoes a sufficient number 
of interactions (> 10) and this makes the generation of UHECRs possible. 




Fig. 4. The UHECR repartition (in % of the total UHECR population) as a function of the UHECR creation time 
which mainly happens between r\, [t = Q on the figures) and about 300 r& (t = 10 4 s) in the Monte-Carlo simulation. 
Top a magnetic field at Tb of 10 4 G : (left At w = 0.5 s, right At w — 10 s). Bottom magnetic field at r\, of 10 5 G (left 
At w — 0.5 s and right At w = 10 s). 



where p+ is reached when the p hotospheric radius is equal 
to ri, ((Gialis fc Pelletier! Eo03). The improvement of the 
acceleration by the process presented in this paper does 
not significantly change the order of magnitude of our pre- 
vious estimate. The number of pp-neutrinos is greater and 
could allow the formation of a spectrum, since we could 
get 10 8 neutrinos under the same conditions. 



6.2.2. Gamma emission by ir -decay 

The pp-process generates 7r°-mesons that decay into 
two 7-photons. This emission deserves more attention be- 
cause of a possible detection in the GeV range, which 
would constitute an interesting diagnosis. For a e~ 2 spec- 
trum of cosmic rays, the gamma flux due to 7r°-decay is 
estimated by 



SV° ~ c a pp n p N± I — 



(21) 
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for e 7 > eo and S^o (e 7 ) = otherwise, where eo — 67 McV 
in the co-moving frame, which corresponds to about 20 
GeV for the observer. Actually the fireball is optically thin 
to Compton scattering at 67 MeV and beyond, thanks to 
the Klein-Nishina regime, and also to gamma-gamma pair 
production for two possible target photons (as this will be 
detailed in a forthcoming paper) . This spectrum thus has a 
power law tail that extends over several decades; the high- 
est part would reveal the presence of UHECRs, but would 
require a very high sensitivity. Detection of gamma-rays 
in that energy range from GRBs seems to be expected by 
ground gamma-ray telescopes, as suggested by the anal- 
ysis of the spectrum of GRB9 70417a by the Milagro col- 
laboration l|Atkins et al 

The 7r°-decay spectrum must be compared with the SSC- 
spectrum. The latter can be written; 



Sssc(t~/) ~ ca T n s N* 



tpeak 



-1/2 



(22) 



for e 7 > esse with e SS c = £ P eak% eak ~ 10 14 eV (for 
e P eak = 1 MeV for the observer and B+{rb) = 10 4 G) the 
energy of the SSC-spectrum is maximum and Sssc( e j) = 
otherwise, where N* is the relativistic fraction of the 
electron population. 

The diagnosis of the 7r°-decay should be fairly easy since 
it radiates at much higher energy than the synchrotron 
emission and at much lower energy than the SSC-emission. 
When rj < 77* /2, the fraction of the energy radiated by ir°- 
decay is quite significant, a few percent; indeed 



E 



■ 0.07 x 



f] 2 t 



10~ 2 - HT 1 . 



(23) 



with to = ro/c ~ 1 ms. 



Such a detection (by HESS observatory, 5@5 exper- 
iment and GLAST) would provide a very nice estimate 
of the baryon load together with evidence of cosmic ray 
generation. 

7. Conclusion 

We have shown through accurate investigation of scat- 
tering properties that the usual Fermi processes are un- 
able to generate UHECRs in GRBs. First of all, it is well 
known that the highly relativistic external shock with a 
standard interstellar medium cannot achieve the goal of 
UHECR generation. The hope that it could do so was 
mostly based on the internal shock model but the irrele- 
vant Bohm's scaling was used; moreover, a slowly decreas- 
ing magnetic field was assumed. Remov ing these extreme 
assum ptions made the goal impossible ()Gialis fc Pelletier] 
l2003h . 

However, because the global energy budget is in favour of 
copious cosmic ray generation, we thought it interesting 
to design an additional Fermi process that would stretch 
the cosmic ray distribution tail. Indeed, the acceleration 
through multi-scattering off magnetized relativistic fronts 



turns out to be more efficient even with a field decreasing 
in r~ 2 . In this paper, we found that a sizeable fraction 
(~ 10~ 5 ) of the cosmic ray population reaches the UHE 
range, which is even more than necessary. 
There is another important physical issue, namely the es- 
cape of the UHECRs from the fireball towards the ob- 
server. The relativistic Fermi process that we proposed 
generates UHECRs in a thin layer located at the very be- 
ginning of the internal shock phase. The question arises 
whether these particles suffer a further expansion loss. In 
fact, because the magnetic field decreases faster than r _1 , 
the interaction energy limit with the magnetic fronts de- 
creases and the UHECRs are no longer scattered by them. 
Thus, because of the lack of scattering, they cannot expe- 
rience the expansion loss, contrarily to models expecting 
B oc r . Therefore, they directly travel across the shells, 
benefiting of the Lorentz boost. For a magnetic field de- 
creasing as r -2 , the integrated spectrum is in e~ 2 over 4 

decades. 

In our previous paper I Gialis fc Pelletier1l2003h . we pre- 
dicted a possible double neutrino emission; even though 
this improvement of the acceleration process does not 
significantly change the estimated fluxes, the collima- 
tion allows detectable events, especially with pp-neutrinos. 
However, we found a possible diagnosis of cosmic ray gen- 
eration in some numerous class of GRBs, having 77 < ??*/2, 
associated with the signature of 7r°-decay in the GeV 7- 
ray spectrum. This seems to be really observable, with a 
few percent of the GRB energy, for GRBs having a suf- 
ficient baryon load. This signature (in the GeV range by 
HESS observatory, 5@5 experiment and GLAST) would 
be easily distinguished from the synchrotron bump and 
the SSC-bump with a flux as high as a few percent of the 
fireball global energy flux. 

This study gives a new chance to GRBs as sources of 
UHECRs. But their contribution to the diffuse back- 
ground spectrum might not be dominant compared to 
the AGN contribution. This is an import ant issue that 
is under intense debate in the literature JSteckedEoOol 
iBerezinskv et al.ll2003UBahcall fc Waxmarjl2003|) . 



Appendix 

From equation [7| we deduce the following expression 

f(e) = Ae-*(5l-Sl), (24) 

where A = n*/(6 ln7 m + T /3o) and Si = Sup{e/7 m , l/2r}, 
6 2 = Inf{e,2r}. 

According to the equation (0) and for a typical magnetic 
field of 10 5 G at r^, we have supposed that j m > 2T. It 
is noteworthy that the following calculation leads to the 
same result for the hypothesis j m < 2r. Let us distinguish 
two cases : the first one for e > j m gives 5\ — e/j m and 
62 = 2r. The distribution function becomes 



8F d - — 



OC e 



(25) 



10 D. Gialis and G. Pelletier: Which acceleration process for UHECRs in GRBs ? 

for e < 2Tj m and / e > 7m (e) = beyond 2T-f m . 

The second case is considered for e < j m and it can be 
divided in two sub cases : first, for e > 2T, the distribution 
function can be simplified and becomes such that 

/2r<e< 7m (e)cxe- 2 . (26) 

Finally, for e < 21", we have S2 — e 3> <5i and the distribu- 
tion function can be written 

/e< 2 r(e)cxe. (27) 
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